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Wc describe the general phenomenon of 'indueed-eharge electro-osmosis' (ICEO) — the 
nonlinear electro-osmotic slip that occurs when an applied field acts on the ionic charge 
it induces around a polarizable surface. Motivated by a simple physical picture, we cal- 
culate ICEO flows around conducting cyHnders in steady (DC), oscillatory (AC), and 
suddenly-applied electric fields. This picture, and these systems, represent perhaps the 
clearest example of nonlinear electrokinetic phenomena. We complement and verify this 
physically-motivated approach using a matched asymptotic expansion to the electroki- 
netic equations in the thin double-layer and low potential limits. ICEO slip velocities 
vary like Ug (x EqL, where Eq is the field strength and i is a geometric length scale, and 
are set up on a time scale Tc = XdL/D, where A_d is the screening length and D is the 
ionic diffusion constant. We propose and analyze ICEO microfiuidic pumps and mixers 
that operate without moving parts under low applied potentials. Similar flows around 
metallic colloids with flxed total charge have been described in the Russian literature 
(largely unnoticed in the West). ICEO flows around conductors with flxed potential, on 
the other hand, have no colloidal analog and offer further possibilities for microfluidic 
applications. 



1. Introduction 

Recent developments in micro-fabrication and the technological promise of microfluidic 
'labs on a chip' have brought a renewed interest to the study of low-Reynolds number 
flows pone fc Kim (2001)| [Whitesides fc Stroock (2001)| [Reyes et al. (2002)| ). The fa- 



miliar techniques used in larger-scale applications for fluid manipulation, which often 
exploit fluid instabilities due to inertial nonlinearities, do not work on the micron scale 
due to the pre-eminence of viscous damping. The microscale mixing of miscible fluids 
must thus occur without the benefit of turbulence, by molecular diffusion alone. For 
extremely small devices, molecular diffusion is relatively rapid; however, in typical mi- 
crofluidic devices with 10-100 ^m features, the mixing time can be prohibitively long (of 
order 100 s for molecules with diffusivity 10~^° m^/s). Another limitation arises because 
the pressure-driven ffow rate through small channels decreases with the third or fourth 
power of channel size. Innovative ideas are thus being considered for pumping, mix- 
ing, manipulating and separating on the micron length scale {e.g. Beebe et al. (2002){ 



Whitesides et al. (2001) I. Naturally, many focus on the use of surface phenomena, owing 
to the large surface to volume ratios of typical microffuidic devices. 

Electrokinetic phenomena provide one of the most popular non-mechanical techniques 
in microffuidics. The basic idea behind electrokinetic phenomena is as follows: locally 
non-neutral fluid occurs adjacent to charged solid surfaces, where a diffuse cloud of 
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oppositely-charged counter-ions 'screen' the surface charge. An externally-apphed electric 
field exerts a force on this charged diffuse layer, which gives rise to a fluid flow relative to 
the particle or surface. Elcctrokinctic flow around stationary surfaces is known as clcctro- 
osniotic flow, and the elcctrokinctic motion of frecly-suspcnded particles is known as elec- 
trophoresis. Electro-osmosis and electrophoresis find wide application in analytical chem- 
istry (|Bruin (2000)|), genomics ( [Landers (2003)| ) and proteomics ( |Figeys fc Pinto (2001)| 



Dolnik fc Hutterer (200T]|. 



The standard picture for electrokinetic phenomena involves surfaces with fixed, con- 
stant charge (or, equivalently, zeta potential C, defined as the potential drop across the 



screening cloud). Recently, variants on this picture have been explored. Anderson (1985) 



demonstrated that interesting and counter-intuitive effects occur with spatially inhomo- 
geneous zeta-potentials, and showed that the electrophoretic mobility of a colloid was 
sensitive to the distribution of surface charge, rather than simply the total net charge. 



Anderson & Idol (1985) explored electro-osmotic flow in inhomogeneously-charged pores, 



and found eddies and recirculation regions. Ajdari (1995; 2001) and Stroock et al. (2000) 



showed that a net electro-osmotic flow could be driven either parallel or perpendicular 
to an applied field by modulating the surface and charge density of a microchannel, and 



Gitlin et al. (2003) have implemented these ideas to make a 'transverse electrokinetic 



pump'. Such transverse pumps have the advantage that a strong field can be estab- 
lished with a low voltage applied across a narrow channel. Long & Ajdari (1998) exam- 



ined electrophoresis of patterned colloids, and found example colloids whose electropho- 
retic motion is always transverse to the applied electric field. Finally, Long et aZ.(1999)1 
Ghosal (2003)1 and others have studied electro-osmosis along inhomogeneously-charged 
channel walls (due to, e.g. , adsorption of analyte molecules), which provides an addi- 
tional source of dispersion that can limit resolution in capillary electrophoresis. 

Other variants involve surface charge densities that are not fixed, but rather are in- 
duced (either actively or passively). For example, the effective zeta potential of chan- 
nel walls can be manipulated using an auxiliary electrode to improve separation ef- 
ficiency in capillary electrophoresis ( Lee et al. (1990)1 Hayes & Ewing (1992) I and, by 
analogy with the electronic field-effect transistor, to set up 'field-effect electro-osmosis' 
( Ghowsi fc Gale (1991)|[Gajar fc Geis (1992)[[Schasfoort et al. (1999) 1. 

Time- varying, inhomogeneous charge double-layers induced around electrodes give rise 
to interesting effects as well. Trau et al. (1997) and Yeh et al. (1997) demonstrated that 



colloidal spheres can spontaneously self-assemble into crystalline aggregates near elec- 
trodes under AC applied fields. They proposed somewhat similar electrohydrodynamic 
mechanisms for this aggregation, in which an inhomogeneous screening cloud is formed by 
(and interacts with) the inhomogeneous applied electric field (perturbed by the sphere), 
resulting in a rectified electro-osmotic flow directed radially in toward the sphere. More 
recently, Nadal et al. (2002a) performed detailed measurements in order to test both the 
attractive (electrohydrodynamic) and repulsive (electrostatic) interactions between the 
spheres, and Ristenpart et al. (2003) explored the rich variety of patterns that form when 
bidisperse colloidal suspensions self-assemble near electrodes. 

A related phenomenon allows steady electro-osmotic flows to be driven using AC elec- 
tric fields. Ramos et al. (1998; 1999) and Gonzalez et al. (2000) theoretically and experi- 
mentally explored 'AC electro-osmosis', in which a pair of adjacent, flat electrodes located 
on a glass slide and subjected to AC driving, gives rise to a steady electro-osmotic flow 
consisting of two counter-rotating rolls. Around the same time, Ajdari (2000) theoreti- 
cally predicted that an asymmetiic array of electrodes with applied AC fields generally 
pumps fluid in the direction of broken symmetry ('AC pumpin g'). [Brown et al. (2001)| 
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Studer et al. (2002)[ and Mpholo et al. (2003)| have since developed AC electrokinetic 



micro-pumps based on this effect, and Ramos et al. (2003) have furthered their analysis. 

Both AC colloidal self-assembly and AC electro-osmosis occur around electrodes whose 
potential is externally controlled. Both effects are strongest when the voltage oscillates 
at a special AC frequency (the inverse of the charging time discussed below), and both 
effects disappear in the DC limit. Furthermore, both vary with the square of the applied 
voltage Vq. This nonlinear dependence can be understood qualitatively as follows: the 
induced charge cloud/zeta potential varies linearly with Vq, and the resulting flow is 
driven by the external field, which also varies with Vq. On the other hand, DC colloidal 
aggregation, as explored by Solomentsev et al. (1997) requires large enough voltages to 



pass a Faradaic current, and is driven by a different, linear mechanism. 

Very recently in microfluidics, a few cases of non-linear electro-osmotic flows around 
isolated and inert (but polarizable) objects have been reported, with both AC and DC 
forcing. In a situation reminiscent of AC electro-osmosis, Nadal et al. (20026) studied the 



micro-flow produced around a dielectric stripe on a planar blocking electrode. In rather 
different situations, Thamida & Chang (2002) observed a DC nonlinear electrokinetic 



jet directed away from a protruding corner in a dielectric microcliannel, far away from 
any electrode, and Takhistov et al. (2003) observed electrokinetically driven vortices near 
channel junctions. These studies (and the present work) suggest that a rich variety of 
nonlinear electrokinetic phenomena at polarizable surfaces remains to be exploited in 
microfluidic devices. 

In colloidal science, nonlinear electro-osmotic flows around polarizable (metallic or di- 
electric) particles were studied almost two decades ago in a series of Ukrainian papers, 
reviewed by Murtsovkin (1996)} that has gone all but unnoticed in the West. Such flows 
occur when the applied field acts on the component of the double-layer charge that has 
been polarized by the field itself. This idea can be traced back at least to Levich (1962)| 
who discussed the dipolar charge double layer (using the Helmholtz model) that is induced 
around a metallic colloidal particle in an external electric field and touched upon the 
quadrupolar flow that would result. Simonov & S. S. Dukhin (1973) calculated the struc- 
ture of the (polarized) dipolar charge cloud in order to obtain the electrophoretic mobility. 



without concentrating on the resulting flow. Gamayunov, Murtsovkin & A. S. Dukhin (1986) 



and A. S. Dukhin & Murtsovkin (1986) first explicitly calculated the nonlinear electro- 
osmotic flow arising from double-layer polarization around a spherical conducting par- 
ticle, and|A. S. Dukhin (1986)] 



extended this calculation to include a dielectric surface 



coating (as a model of a dead biological cell). Experimentally, Gamayunov et al. (1992) 



observed a non-linear flow around spherical metallic colloids, albeit in a direction oppo- 
site to predicitions for all but the smallest particles. They argued that a Faradaic current 
(breakdown of ideal polarizibility) was responsible for the observed flow reversal. 

This work followed naturally from many earlier studies on 'non-equilibrium electric 
surface phenomena' reviewed by S. S. Dukhin (1993)[ especially those focusing on the 'in- 



duced dipole moment' (IDM) of a colloidal particle reviewed by S. S. Dukhin & Shilov (1980) 
Following Overbeek (1943)| who was perhaps the first to consider non-uniform polariza- 
tion of the double layer in the context of electrophoresis, S. S. Dukhin (1965)] S. S. Dukhin fc Semenikhin (1970)[ 
andS. S. Dukhin & Shilov (1974) predicted the electrophoretic mobility of a highly charged 
non-polarizable sphere in the thin double-layer limit, in good agreement with the later 
numerical solutions of ( 



O'Brien fc White (1978) (See, e.g., Lyklema (1991) ) In that case. 



diffuse charge is redistributed by surface conduction, which produces an IDM aligned with 
the field and some variations in neutral bulk concentration, and secondary electro-osmotic 
and diffusio-osmotic flows develop as a result. Shilov & S. S. Dukhin (1970) extended this 



work to a non-polarizable sphere in an AC electric field. Simonov & S. S. Dukhin (1973) 
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and Simonov & Shilov (1973) performed similar calculations for an ideally polarizable, 

conducting sphere, which typically exhibits an IDM opposite to the applied field. [Simonov fc Shilov (1977) 



revisited this problem in the context of dielectric dispersion and proposed a much simpler 



RC-circuit model to explain the sign and frequency dependence of the IDM. O'Brien & White (1978) 



performed a numerical solution of the full ion, electrostatic, and fluid equations with ar- 
bitrarily thick double-layer, which naturally incorporated the effects of double-layer po- 
larization. O'Brien & Hunter (1981) and O'Brien (1983) following Dukhin's approach, 
arrived at a simpler, approximate expression that incorporated double-layer polariza- 
tion in the thin double-layer limit, and compared favorably with the numerical calcula- 
tions of O'Brien & White (1978) Nonetheless, it seems a detailed study of the associ- 
ated electro-osmotic flows around polarizable spheres did not appear until the paper of 



Gamayunov, Murtsovkin & A. S. Dukhin (1986) cited above. 



In summary, we note that electrokinetic phenomena at polarizable surfaces share a fun- 
damental feature: all involve a nonlinear flow component in which double-layer charge 
induced by the applied field is driven by that same field. To emphasize this common 
mechanism, we suggest the term 'induced-charge electro-osmosis' (ICEO) for their de- 
scription. Specific realizations of ICEO include AC electro-osmosis at micro-electrodes, 
AC pumping by asymmetric electrode arrays, DC electrokinetic jets around dielectric 
structures, DC and AC flows around polarizable colloidal particles, and the situations 
described below. Of course, other electrokinetic effects may also occur in addition to 
ICEO in any given system, such as those related to bulk concentration gradients pro- 
duced by surface conduction or Faradaic reactions, but we ignore such complications here 
to highlight the basic effect of ICEO. 

In the present work, we build upon this foundation of induced-charge electrokinetic 
phenomena, with a speciflc eye towards microfluidic applications. ICEO flows around 
metallic colloids, which have attraced little attention compared to non-polarizable ob- 
jects of fixed zeta potential, naturally lend themselves for use in microfluidic devices. 
In that setting, the particle is replaced by a fixed polarizable object which pumps the 
fluid in response to applied fields, and a host of new possibilities arise. The ability to 
directly control the position, shape, and potential of one or more 'inducing surfaces' in 
a microchannel allows for a rich variety of effects that do not occur in colloidal systems. 
Before we begin, we note the difference between ICEO and 'electrokinetic phenomena of 



the second kind', reviewed by S. S. Dukhin (1991) and studied recently by Ben & Chang (2002) 



in the context of microfiuidic applications. Significantly, second- kind electrokinetic effects 
do not arise from the double-layer, being instead driven by space charge in the bulk solu- 
tion. They typically occur around ion-selective porous granules subject to applied fields 
large enough to generate strong currents of certain ions through the liquid/solid interface. 
This leads to large concentration variations and space charge in the bulk electrolyte on one 
side, which interact with the applied field to cause motion. Barany et al. (1998) studied 
the analogous effect for non-porous metallic colloids undergoing electrochemical reac- 
tions at very large Faradaic currents (exceeding diffusion limitation). In contrast, ICEO 
occurs around inert polarizable surfaces carrying no Faradaic current in contact with a 
homogeneous, quasi-neutral electrolyte and relies on relatively small induced double-layer 
charge, rather than bulk space charge. 

The article is organized as follows: ^provides a basic background on electrokinetic 
effects, and 33 develops a basic physical picture of induced-charge electro-osmosis via 
calculations of steady ICEO fiow around a conducting cylinder. SQl examines the time- 
dependent ICEO fiow for background electric fields which are suddenly-applied f ^4.1|l or 
sinusoidal ft l4.2(l . SjHl describes some basic issues for ICEO in microfluidic devices, such 
as coupling to the external circuit f tl5.1|l and the phenomenon of fixed-potential ICEO 
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Figure 1. Left: A charged solid surface in an electrolytic solution attracts an oppositely-charged 
'screening cloud' of width ~ Ad. An electric field applied tangent to the charged solid surface 
gives rise to an electro-osmotic flow 12.41 . Right: An electric field applied to an electrolytic 
solution containing a suspended solid particle gives rise to particle motion called electrophoresis, 
with velocity equal in magnitude and opposite in direction to 12.411 . 



( ^5.211 . Some specific designs for microfluidic pumps, junction switches, and mixers are 
discussed and analyzed in t |5.3l ^ investigates the detrimental effect of a thin dielectric 
coating on a conducting surface and calculates the fCEO flow around non-conducting 
dielectric cylinders. Ogives a systematic derivation of fCEO in the limit of thin double 
layers and small potentials, starting with the basic electrokinetic equations and employ- 
ing matched asymptotic expansions, concluding with a set of effective equations (with 
approximations and errors quantified) for the fCEO flow around an arbitrarily-shaped 
particle in an arbitrary space- and time-dependent electric field. The interesting conse- 
quences of shape and field asymmetries, which generally lead to electro-osmotic pumping 
or electrophoretic motion in AC fields, are left for a companion paper. The reader is 
referred to Bazant & Squires (2004) for an overview of our results. 



2. Classical ('fixed- charge') electro-osmosis 

Electrokinetic techniques provide some of the most popular small-scale non-mechanical 
strategies for manipulating particles and fluids. We present here a very brief introduction. 
More detailed accounts are given by Lyklema (1991),|Ifunter (2000) and Russel et al. (1989) 



2.1. Small zeta potentials 

A surface with charge density q in an aqueous solutions attracts a screening cloud of 
oppositely-charged counter- ions to form the electrochemical 'double layer', which is ef- 
fectively a surface capacitor. In the Debye-Hiickel limit of small surface charge, the excess 
diffuse ionic charge exponentially screens the electric field set up by the surface charge 
(figure^), giving an electrostatic potential 

(/) = -^e-'"^ = Ce-''^ (2.1) 

Here £w ~ SOeq is the dielectric permittivity of the solvent (typically water) and Eq is 
the vacuum permittivity. The 'zeta potential', defined by 
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reflects the electrostatic potential drop across the screening cloud, and the Debye 'screen- 
ing length' is defined for a symmetric z:z electrolyte by 



with bulk ion concentration tioi (monovalent) ion charge e, Boltzmann constant fcs and 
temperature T. Because the ions in the diffuse part of the double layer are approximately 
in thermal equilibrium, the condition for a 'small' charge density (or zcta potential) is 
C,-€.kBT/ze. 

An externally applied electric field exerts a body force on the electrically charged 
fluid in this screening cloud, driving the ions and the fluid into motion. The resulting 
electro-osmotic fluid flow (figure QJi) appears to 'slip' just outside the screening layer 
of width Xd- Under a wide range of conditions, the local slip velocity is given by the 
Helmholtz-Smoluchowski formula, 

u, = -^E||, (2.4) 
r] 

where rj is the fluid viscosity and E|| is the tangential component of the bulk electric 
field. 

This basic electrokinetic phenomenon gives rise to two related effects, electro-osmosis 
and electrophoresis, both of which find wide application in analytical chemistry, microfiu- 
idics, colloidal self-assembly, and other emerging technologies. Electro-osmotic flow occurs 
when an electric field is applied along a channel with charged walls, wherein the electro- 
osmotic slip at the channel walls gives rise to plug fiow in the channel, with velocity 
given by H2.4|l . Because the electro-osmotic flow velocity is independent of channel size, 
(in contrast to pressure-driven fiow, which depends strongly upon channel size), electro- 
osmotic pumping presents a natural and popular technique for fluid manipulation in small 
channels. On the other hand, when the solid/fluid interface is that of a freely-suspended 
particle, the electro-osmotic slip velocity gives rise to motion of the particle itself (figure 
^), termed electrophoresis. In the thin double-layer limit, the electrophoretic velocity is 
given by Smoluchowski's formula, 

U=^Eoo = MeEoo, (2.5) 

V 

where Eqo is the externally-applied field, and /ie — s^(^/ri is the electrophoretic mobility 
of the particle. 

2.2. Large zeta potentials 

For 'large' zeta potentials, comparable to or exceeding ksT/ze, the exponential profile 
of the diffuse charge 1)2. 1() and the linear charge- voltage relation 1)2. 2() are no longer valid, 
but the diffuse part of the double layer remains in thermal equilibrium. As a result, the 
potential in the diffuse layer satisfles the Poisson-Boltzmann equation. For symmetric, 
binary electrolyte, the classical Gouy-Chapman solution yields a nonlinear charge- voltage 
relation for the double layer, 

q{C) = AnnzeXo sinh (^^^) , (2-6) 

in the thin double-layer limit. This relation may be modified to account for microscopic 
surface properties, such as a compact Stern layer, a thin dielectric coating, or Faradaic 
reactions, which enter via the boundary conditions to the Poisson-Boltzmann equation. 
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The key point here is simply that q generaUy grows exponentiaUy with ze^/ksT, which 
has important imphcations for time-dependent problems involving double-layer relax- 
ation, as reviewed by Bazant, Thornton fc Ajdari (2004)1 

Remarkably, the Hclmholtz-Smoluchowski tormuia l|2.4|l for the electro-osmotic slip 
remains valid in the nonlinear regime, as long as 



A 



D 



exp 



ze(^ 
2kBT 



« 1, 



(2.7) 



where a is the radius of curvature of the surface (Hunter 2000). When (|2.7|) is violated, 
ionic concentrations in the diffuse layer differ significantly from their bulk values, and 
surface conduction through the diffuse layer becomes important. As a result, the electro- 
phoretic mobility, /ig, becomes a nonlinear function of the dimensionless ratio 

Du=— , (2.8) 



of the surface conductivity, CTs, to the bulk conductivity, a. Though this ratio was first 
noted by Bikerman (1940)| we follow Lyklema (1991) in referring to H2.8|l as the 'Dukhin 
number', in honor of Dukhin's pioneering calculations of its effect on electrophoresis. We 
note also that essentially the same number has been called a, (3, and A by various authors 



in the Western literature, and 'Rel' by S. S. Dukhin (1993) 



Bikerman (1933; 1935) made the first theoretical predictions of surface conductance, 
(Ts, taking into account both electromigration and electro-osmosis in the diffuse layer. 
The relative importance of the latter is determined by another dimensionless number 

(2.9) 



ze 



2e^ 
rjD 



Using the result of Deryagin fc S. S. Dukhin (1969)| Bikerman's formula for takes a 
simple form for a symmetric binary electrolyte, yielding 



Du = 



2Ac(l + m) 



cosh 



zeC 



- 1 



= 4Mi±Z!^sinh2 



zeC, 
^keT 



(2.10) 



for the Dukhin number. In the limit of infinitely thin double layers, where (j2.7|l holds, 
the Dukhin number vanishes, and the electrophoretic mobility tends to Smoluchowski's 
value H2.4() . For any finite double-layer thickness, however, highly charged particles (with 
C, > ksT/ze) are generally described by a non-negligible Dukhin number, and surface 
conduction becomes important. This leads to bulk concentration gradients and a non- 
uniform diffuse-charge distribution around the particle in an applied field, which modifies 
its electrophoretic mobility, via diffusiophoresis and concentration polarization. For more 



details, the reader is referred to Lyklema (1991) and S. S. Dukhin (1993) 



3. Induced-charge electro-osmosis: fundamental picture 

The standard electrokinetic picture described above involves the interaction of an ap- 
plied field and a surface of fixed charge, wherein the electro-osmotic fiow is linear in 
the applied field. Here we focus on the nonlinear phenomenon of ICEO at a polarizable 
(metal or dielectric) surface. As a consequence of nonlinearity, ICEO can be used to 
drive steady electro-osmotic flows using AC or DC fields. The nonlinearity also allows for 
larger fluid velocities and a richer, geometry-dependent flow structure. These properties 
stand in stark contrast to classical electro-osmosis described above, which, e.g. , gives 
zero time-averaged fiow in an AC field. 
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Figure 2. The evolution of the electric field around a solid, ideally polarizable conducting 
cylinder immersed in a liquid electrolyte, following the imposition of a background DC field 
at t = G (a), where the field lines intersect normal to the conducting surface. Over a charging 
time Tc = Xoa/D, a dipolar charge cloud forms in response to currents from the bulk, reaching 
steady state (b) when the bulk field profile is that of an insulator. The resulting zeta potential, 
however, is nonuniform. 



This section gives a physically clear (as well as quantitatively accurate) sense of ICEO 
in perhaps the simplest case: an ideally polarizable metal cylinder in a suddenly ap- 
plied, uniform electric field. This builds on the descriptions of double-layer polariza- 
tion for an uncharged metallic sphere by Levich (1962)1 [Simonov fc Shilov (1973)1 and 

Simonov & Shilov (1977) and the associated steady ICEO flow by Gamayunov, Murtsovkin & A. S. Dukhin (1986) 



We postpone for §7| a more detailed and general analysis, justifying the approxima- 
tions made here for the case of thin double layers (Ad <C a) and weak applied fields 
(E'ofl ksT/ze). Since Du <C 1 in these limits, surface conduction and concentration 
polarization can be safely ignored, and ICEO becomes the dominant clcctrokinctic effect 
around any inert, highly polarizable object. 



3.1. Steady ICEO around an uncharged conducting cylinder 

The basic phenomenon of ICEO can be understood from figures [2 and |2| Immediately 
after an external field E = EqZ is applied, an electric field is set up so that field lines 
intersect conducting surfaces at right angles (figure Et) . Although this represents the 
steady-state vacuum field configuration, mobile ions in electrolytic solutions move in 
response to applied fields. A current J ~ ctE drives positive ions into a charge cloud 
on one side of the conductor (z < 0), and negative ions to the other (z > 0), inducing 
an equal and opposite surface charge on the conducting surface. A dipolar charge cloud 
grows so long as a normal field injects ions into the induced double layer, and steady 
state is achieved when no field lines penetrate the double-layer (figurc|2h) ■ The tangential 
field i?| [ drives an electro-osmotic slip velocity (|2.4ll proportional to the local double-layer 
charge density, driving fluid from the 'poles' of the particle towards the 'equator' (figure 
I2K)- An AC field drives an identical flow, since an oppositely-directed field induces an 
oppositely-charged screening cloud, giving the same net flow. The ICEO flow around a 
conducting cylinder with non-zero total charge, shown in figurc|2t), simply superimposes 
the nonlinear ICEO flow of figure with the usual linear electro-osmotic flow. 
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Figure 3. The steady-state induced-charge electro-osmotic flow around (a) a conducting cylin- 
der with zero net charge and (b) a positively charged conducting cylinder. The ICEO slip velocity 
depends on the product of the steady field and the induced zeta potential. The flow around an 
uncharged conducting cylinder (a) can thus be understood qualitatively from figure|2h, whereas 
the charged sylinder (b) simply involves the superposition of the standard electro-osmotic flow. 



As a concrete example for quantitative analysis, we consider an isolated, uncharged 
conducting cylinder of radius a submerged in an electrolyte solution with very small 
screening length A_d ^ a. An external electric field E^^i is suddenly applied at < = 0, and 
the conducting surface forms an equipotential surface, giving a potential 



a2 



E^zil-^]. (3.1) 



r 

Electric field lines intersect the conducting surface at right angles, as shown in figurcEt. 

Due to the electrolyte's conductivity a, a non-zero current J = crE drives ions to the 
cylinder surface. In the absence of electrochemical reactions at the conductor/electrolyte 
interface {i.e. at sufRciently low potentials that the cylinder is 'ideally polarizable'), mo- 
bile solute ions accumulate in a screening cloud adjacent to the solid/liquid surface, 
attracting equal and opposite 'image charges' within the conductor itself. Thus the con- 
ductor's surface charge density q — induced by the growing screening cloud — changes 
in a time-dependent fashion, via 

«.j.f^.E.f. (3.2) 

Using the linear relationship (|2.2|l between surface charge density and zeta potential, this 
can be expressed as 

^ = ^E.r. (3.3) 

A dipolar charge cloud grows, since positively-charged ions are driven into the charge 
cloud on the side of the conductor nearest the field source (z < in this case), and 
negatively-charged ions are driven into the charge cloud on the opposite side. As ions are 
driven into the screening charge cloud, field lines are expelled and the ionic flux into the 
charge cloud is reduced. 

The system reaches a steady state configuration when all field lines are expelled (f • 
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E(a) = 0). This occurs when the electrostatic potential outside of the charge cloud is 
given by 

shown in figure \3p. The steady-state electrostatic configuration is thus equivalent to 
the no-flux electrostatic boundary condition assumed in the analysis of 'standard' elec- 
trophoresis. In the present case, however, the steady-state configuration corresponds to 
a cylinder whose zeta potential varies with position according to 

C^{e) ^ (Po- <f)f {a) :^2EQa cos 6, (3.5) 

where we assume the conductor potential (po to vanish. While the steady-state electric 
field has no component normal to the charge cloud, its tangential component, 

^•E = -2£;osin6l, (3.6) 

drives an induced-charge electro-osmotic flow, with slip velocity given by (E3I). Now, 
however, the (spatially- varying) surface potential C, is given by (|3.5|) . Because the charge 
cloud is itself dipolar, the tangential field drives the two sides of the charge cloud in 
opposite directions — each side away from the poles — resulting in a quadrupolar electro- 
osmotic slip velocity 

u, = 2[/osin26l ^, (3.7) 
where C/q is the natural velocity scale for ICEO, 

f/o = (3.8) 

One power of sets up the 'induced-charge' screening cloud, and the second drives the 
resultant electro-osmotic flow. 

The fluid motion in this problem is reminiscent of that around a fluid drop of one 
conductivity immersed in a fluid of another conductivity subjected to an external electric 



field, studied by Taylor (1966) By analogy, we find the radial and azimuthal fluid velocity 



components of the fluid flow outside of the cylinder to be 

/ 2 2 \ 3 
Ur = 2 "^" 7^ V ocos26>, 2^[/osin26l. (3.9) 

For comparison, analogous results for the steady-state ICEO flow around a sphere, some 
of which were derived by Gamayunov, Murtsovkin fc A. S. Dukhin (1986)| are given in 



Tabled 

Although we focus on the limit of linear screening in this paper, (|3.9() should also hold 
for non-linear screening (^ w ksT/ze) in the limit of thin double layers, as long as the 
Dukhin number remains small and (|2.7|l is satisfied. The relevant zeta potential in these 
conditions, however, is not the equilibrium zeta potential (Co = 0) but the typical induced 
zeta potential, C sa Eija, which is roughly the applied voltage across the particle. 

Finally, although we have specifically considered a conducting cylinder, a similar pic- 
ture clearly holds for more general shapes. More generally, ICEO slip velocities around 
arbitrarily-shaped inert objects in uniform applied fields are directed along the surface 
from the 'poles' of the object (as defined by the applied field), towards the object's 
'equator'. 

3.2. Steady ICEO around a charged conducting cylinder 

Until now, we have assumed the cylinder to have zero net charge for simplicity. A cylinder 
with non-zero equilibrium charge density qo = Q/Ana^, or zeta potential, Co = EuiKqo, in 
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Cylinder 



Sphere 



Initial potential (j>i 

Steady-state potential <ps 
Steady-state zeta potential 
Radial flow Ur 
Azimuthal flow ug 
Charging Timescale 

Induced dipole strength: 

g{t) for suddenly applied field Eo 

g for AC field Re(Eoe""'*) 



Eoz 1 - ^ 



2Eoacos9 



if/ocos 26 



^(7o sin 261 

_ Xpa 

g{t) = 1 - 2e-*/-= 



l-\-iuiTc 



~Eoz (^1 - ^ 
-Eoz + 4 



lEoa cos t 



Uo (1 + 3cos 261) 



rf/osin 261 



ffW = i(l-3e-2*/-) 



Table 1. Electrostatic and hydrodynamic quantities for the induced-charge electro-osmotic 
(ICEO) fiow around conducting spheres and cylinders, each of radius a. Here Uo = ewE pa/ri 
is a characteristic velocity scale, and the induc ed dipole strength g is defined in (14.111 . See 
[Gamayunov, Murtsovkin fc A. S. Dukhin (1986)| for flows around metal colloidal spheres in 
steady and AC fields. 



a suddenly applied field approaches a steady-state zeta-potential distribution, 

C{0) ^ Ca + '^Eoacose, (3.10) 

which has the same induced component in (|3.5|l added to the constant equilibrium value. 
This follows from the linearity of (|3.3|) with the initial condition, C(0,t = 0) = Co- As 
a result, the steady-state electro-osmotic slip is simply a superposition of the 'standard' 
electro-osmotic flow due to the equilibrium zeta potential Coi 

uQ = u3-2^sin0^, (3.11) 
V 

where Ug is the ICEO slip velocity, given in H3.7|l . The associated Stokes flow, which is a 
superposition of the ICEO flow and the 'standard' electro-osmotic flow, and is shown in 
figure I^Jb) . 

The electrophoretic velocity of a charged conducting cylinder can be found using the 
results of [Stone fc Samuel (1996) from which it follows that the velocity of a cylinder 



with prescribed slip velocity Ug{9), but no externally-applied force, is given by the surface- 
averaged velocity, 

1 f^^ 

U= -— / u,{e)dO. (3.12) 
271- Jo 
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The ICEO component (|3.7|l has zero surface average, leaving only the 'standard' electro- 
osmotic slip velocity This was pointed ou t by|Levich (1962)| using a Helm holtz 



model for the induced double layer, and later by Simonov fc S. S. Dukhin (1973)| usmg 
a double-layer structure found by solving the electrokinetic equations. The conducting 
cylinder thus has the same electrophorctic mobility = ^wCo/v as an object of fixed 
uniform charge density and constant zeta potential. This extreme case illustrates the 
result of O'Brien & White (1978) that the electrophorctic mobility does not depend on 
electrostatic boundary conditions, even though the flow around the particle clearly does. 

As above, the steady-state analysis of ICEO for a charged conductor is unaffected by 
nonlinear screening, as long as (|2.7|) is satisfied (and Du <C 1), where the relevant ( is the 
maximum value, C = |Co| + |^oa|, including both equilibrium and induced components. 



4. Time-dependent ICEO 

A significant feature of ICEO flow is its dependence on the square of the electric 
field amplitude. This has important consequences for AC fields: if the direction of the 
electric field in the above picture is reversed, so are the signs of the induced surface 
charge and screening cloud. The resultant ICEO flow, however, remains unchanged: the 
net flow generically occurs away from the poles, and towards the equator. Therefore, 
induced-charge electro-osmotic flows persist even in an AC applied fields, so long as the 
frequency is sufficiently low that the induced-charge screening clouds have time to form. 

AC forcing is desirable in microfluidic devices, so it is important to examine the time- 
dependence of ICEO flows. As above, we explicitly consider a conducting (ideally polar- 
izable) cylinder and simply quote the analogous results for a conducting sphere in Table 
^ Although we perform calculations for the more tractable case of linear screening, we 
briefly indicate how the analysis would change for large induced zeta potentials. Two sit- 
uations of interest are presented: the time-dependent response of a conducting cylinder 
to a suddenly-applied electric held f ti4.1|) and a sinusoidal AC electric field f H4.2|l . We 
also comment on the basic time scale for ICEO flows. 

4.1. ICEO around a conducting cylinder in a suddenly- applied DC field 

Consider first the time-dependent response of an uncharged conducting cylinder in an 
electrolyte when a uniform electric field E = EqZ is suddenly turned on at < = 0. The 
dipolar nature of the external driving suggests a bulk electric field of the form 

cj){r,t) = -Eoz(^l + g{t)^y (4.1) 

so that initially g(0) = —1 H3.1|l . and in steady state g{t — > oo) = 1 H3.4|l . The potential 
of the conducting surface itself remains zero, so that the potential drop across the double 
layer is given by 

cj>{a,9,t) = -C{9) = -^. (4.2) 

Here, as before, we take q to represent the induced surface charge, so that the total charge 
per unit area in the charge cloud is —q. The electric field normal to the surface, found 
from (|4.1(l . drives an ionic current 

J^ = -m^-'^EoCose{l~g), (4.3) 

into the charge cloud, locally injecting a surface charge density q per unit time. We 
express the induced charge density q in terms of the induced dipole g by substituting 
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l|4.1(l into (|4.2(l . take a time derivative, and equate the result with q given by 1)4. 3|1 . This 
results in an ordinary differential equation for the dipole strength g, 

whose solution is 

5(^) = l-2e-*/^^ (4.5) 
Here Tc is the characteristic time for the formation of induced-charge screening clouds, 



D 



(4.6) 



where the definitions of conductivity [a ~ 2noe^D/kBT) and screening length 12.3|l have 
been used. 

The induced-charge screening cloud (with equivalent zeta potential given by H4.2|l ) is 
driven by the tangential field (derived from H4.1|l 'l in the standard way H2.4|) . resulting in 
an induced-charge electro-osmotic slip velocity 

u, = 2C/osin26l (^1 - e"*/^"^)^^. (4.7) 

More generally, the time-dependent slip velocity around a cylinder with a nonzero fixed 
charge (or equilibrium zeta potential (^q) can be found in similar fashion, and results in 
the standard ICEO slip velocity (|4.7|l with an additional term representing standard 
electro-osmotic slip (|3.11|) 

uf ^Us-2^sme(l-e-'/-A9. (4.8) 



V 

Note that 1)4.8(1 grows more quickly than 1)4.7(1 initially, but that ICEO slip eventually 
dominates in strong fields, since it varies with Eq, versus Eq for the standard electro- 
osmotic slip. 

4.2. ICEO around a conducting cylinder in a sinusoidal AC field 

An analogous calculation can be performed in a straightforward fashion for sinusoidal 
applied fields. Representing the electric field using complex notation, E = Eqc^'^^z, where 
the real part is implied, we obtain a time-dependent zeta potential 

C = 2£;oa cos 61 Rc , (4.9) 

V 1 + tUJTc ) 

giving an induced-charge electro-osmotic slip velocity 



Us = 21}^ sin 2B 
with time-averaged slip velocity 



Rc 



1 + iUTc 



Uosm29 



(4.10) 



(4.11) 



In the low- frequency limit luTc 1, the double-layer fully develops in phase with the 
applied field. In the high-frequency limit luTc ^ 1, the double-layer does not have time 
to charge up, attaining a maximum magnitude (ljTc)^^ with a tt/2 phase shift. Note 
that this analysis assumes that the double-layer changes quasi-steadily, which requires 
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4.3. Time scales for IC'EO flows 

Before continuing, it is worth emphasizing the fundamental time scales arising in ICEO. 
The basic charging time exceeds the Debye time for diffusion across the double-layer 
thickness, = ^jj/D = Swja^ by a geometry-dependent factor, a/Xo, that is typically 
very large, is also much smaller than the diffusion time across the particle, = a} jD. 
The appearance of this time scale for induced dipole moment of a metallic colloidal 



sphere has been explained by Simonov & Shilov (1977) using a simple RC-circuit analogy, 
consistent with the detailed analysis of Simonov & Shilov (1973) The same time scale, 
Tc, also arises as the inverse frequency of AC electro-osmosis | |Ramos et al. (1999~ l or 
AC pumping ( Ajdari (2000) I at a micro-electrode array of characteristic length a, where 
again an RC-circuit analogy has been invoked to explain the charging process. This simple 
physical picture has been criticized by Scott et al. (2001)1 but Ramos et al. (2001) and 
Gonzalez et al. (2000) have convincingly defended its validity, as in the earlier Russian 



papers on polarizable colloids. 

Although it is apparently not well known in microfluidics and colloidal science, the 
time scale for double-layer relaxation was debated and analyzed in electrochemistry 
in the middle of the last century, after Ferry (1948) predicted that r/) should be the 



charging-time for the double layer at an electrode in a semi-infinite electrochemical cell. 
Buck (1969) explicitly corrected Ferry's analysis to account for bulk conduction, which 
introduces the macroscopic electrode separation a. The issue was definitively settled by 



Macdonald (1970) who explained the correct charging time, Tc, as the 'i?C time' for the 
double-layer capacitor, C = Ew/Xd, coupled to a bulk resistor, R = a/a. Similiar ideas 
were also developed independently a decade later by Kornyshev & Vorotyntsev (1981) in 
the context of solid electrolytes. 

Ferry's model problem of a suddenly imposed surface charge-density in a semi-infinite 
electrolyte (as opposed to a suddenly imposed voltage or background field in a finite 
system) persists in recent textbooks on colloidal science, such as Hunter (2000) and 
Lyklema (1991)[ and only the time scales tq and are presented as relevant for double- 
layer relaxation. This is quite reasonable for non-polarizablc colloidal particles, but we 
stress that the intermediate RC time scale, = ^tdTq, plays a central role for polarizable 
objects that exhibit significant double-layer charging. 

We also mention nonlinear screening effects at large applied fields or large total charges, 
where the maximum total zeta potential, « C,q + E^a, exceeds ksT / ze. The analysis of 
this section can be generalized to account for the 'weakly nonlinear' limit of thin double 
layers, where C, > ksT/ze, but (|2.7|l is still satisfied (and thus Du <C 1 as well). In the 
absence of surface conduction, (|3.2|l still describes double-layer relaxation, but a nonlinear 
charge- voltage relation, such as H2.6|l from Gouy-Chapman theory, must be used. In that 
case, the time-dependent boundary condition (|3.3|) on the cylinder is replaced by 



crE ■ 



where 



Cd(C) 



A, 



cosh 



zeC. 
2kBT 



(4.12) 



(4.13) 



is the nonlinear differential capacitance of the double layer. 

If the induced component of the zeta potential is large, due to a strong applied field, 
E^a > ksT/ze, the charging dynamics in (|4.12() are no longer analytically tractable. Since 
the differential capacitance Cd(C) increases with |C| in any thin double-layer model, the 
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'poles' of the cylinder along the applied field charge more slowly than the sides. However, 
the steady-state field is the same as for linear screening, so long as Du <C 1 . 

If the applied field is weak, but the total charge is large (Co > ksT/ze), (|4.12|) may 
be linearized to obtain the same polarization dynamics and ICEO flow as for Du <C 1. 
However, the 'i?C time', 

T,(Co) =RCd = t, cosh , (4.14) 



increases nonlinearly with (^q, as shown by Simonov & Shilov (1977) The same time 
constant, with Co replaced by Co + Eoa, also describes the late stages of relaxation in 
response to a strong applied field. 

Some implications for ICEO at large voltages in the 'strongly nonlinear' limit of thin 
double layers, where condition (|2.7() is violated and Du ^ 1, are discussed in ^7.5\ In 
this regime, the simple circuit approximation breaks down due to bulk diffusion, and 
secondary relaxation occurs at the slow time scale, Tq = / D. For an interdisciplinary 
review and detailed analysis of double-layer relaxation (without surface conduction or 



flow), the reader is referred to Bazant, Thornton & Ajdari (2004) 



Finally, we note that the oscillatory component of ICEO flows may not obey the quasi- 
steady Stokes equations, due to the finite time scale, r„ = a} jv, for the diffusion of fluid 
vorticity. (Here v = rj/ pis the kinematic viscosity and p the fluid density.) It is customary 
in microfluidic and colloidal systems to neglect the unsteady term, pdu/dt, in the Stokes 
equations, because ions diffuse more slowly than vorticity by a factor of D/v « 10"'^. 
However, the natural time scale for the AC component of AC electro-osmotic flows is 
Tc = ^dci/D, so the importance of the unsteady term in the Stokes equations is governed 
by the dimensionless parameter, 

^ = -^^. (4.15) 

This becomes significant for sufficiently thin double layers, Au/a « 10"'^. Therefore, 
in AC electro-osmosis and other ICEO phenomena with AC forcing at the charging 
frequency, Wc = t'c^^, vorticity diffusion confines the oscillating component of ICEO flow 
to an oscillatory boundary layer of size ^ vX^aj D. However, the steady component of 
ICEO flows is usually the most important, and obeys the steady Stokes equations. 



5. ICEO in microfluidic devices 

We have thus far considered isolated conductors in background fields applied 'at in- 
finity', as is standard in the colloidal context. The further richness of ICEO phenomena 
becomes apparent in the context of microfluidic devices. In this section, we consider 
ICEO in which the external field is applied by electrodes with finite, rather than infinite, 
separation. Furthermore, microfluidic devices allow additional techniques not available 
for colloids: the 'inducing surface' can be held in place and its potential can be actively 
controlled. This gives rise to 'fixed-potential' ICEO, which is to be contrasted with the 
'fixed-total-charge' ICEO studied above. Finally, we present a series of simple ICEO- 
based microfluidic devices that operate without moving parts in AC flelds. 

5.1. Double-layer relaxation at electrodes 

As emphasized above, one must drive a current Jo(0 — <^Eo{t) to apply an electric 
fleld Eo{t) in an homogeneous electrolyte. The electrochemical reactions associated with 
steady Faradaic currents may cause fluid contamination by reaction products or elec- 
trodeposits, unwanted concentration polarization, or permanent dissolution (and thus 
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irreversible failure) of microelectrodes. Therefore, oscillating voltages and non-Faradaic 
displacement currents at 'blocking' electrodes are preferable in microfluidic devices. In 
this case, however, one must take care that diffuse-layer charging at the electrodes does 
not screen the field. 

We examine the simplest case here, which involves a device consisting of a thin con- 
ducting cylinder of radius a ^ \d placed between flat, inert electrodes separated by 
2L ^ 2a. The cylinder is electrically isolated from the rest of the system, so that its 
total charge is fixed. Under a suddenly applied DC voltage, 2Vo, the bulk electric field, 
Eo{t), decays to zero as screening clouds develop at the electrodes. For weak potentials 
(Vq ^ ksT/e) and thin double-layers (A^ a L), the bulk field decays exponentially 

Eo{t) = ^e''/^', (5.1) 



with a characteristic electrode charging time 



D 



(5.2) 



analogous to the cylinder's charging time (|4.6|) . This time-dependent field i?o(i) then acts 
as the 'applied field at infinity' in the ICEO slip formula, H4.7|l . 

The ICEO flow around the cylinder is set into motion exponentially over the cylinder 
charging time, ~ Xua/D, but is terminated exponentially over the (longer) electrode 
charging time, = \dL/D as the bulk field is screened at the electrodes. This interplay 
between two time scales — one set by the geometry of the inducing surface and another 
set by the electrode geometry — is a common feature of ICEO in microfluidic devices. 

This is clearly seen in the important case of AC forcing by a voltage, T^) cos(Li;t), in 
which the bulk electric fleld is given by 

E,(t) = ^Rc ( _^^,-^^t\ (5 3) 



L V 1 + iijJTf. 

Electric fields persist in the bulk solution when the driving frequency is high enough 
(wTg ^ 1) that induced double-layers do not have time to develop near the electrodes. 
Induced-charge electro-osmotic flows driven by applied AC fields can thus persist only 
in a certain band of driving frequencies, < w < r~^, unless Faradaic reactions 
occur at the electrodes to maintain the bulk fleld. In AC electro-osmosis at adjacent sur- 



face electrodes (Ramos et al. (1999) I or AC pumping at an asymmetric electrode array 



( Ajdari (2000) I, the inducing surfaces are the electrodes, and so the two time scales coin- 



cide to yield a single characteristic frequency uJc = l/rg. (Note that Ramos et al. (1999) 



and Gonzalez et al. (2000) use the equivalent form uj ~ aXo/^wL.) 

Table 121 presents typical values for ICEO flow velocities and charging time scales for 
some reasonable microfluidic parameters. For example, an applied electric fleld of strength 
100 V/cm across an electrolyte containing a 10 fj,m. cylindrical post gives rise to an ICEO 
slip velocity of order 1 mm/s, with charging times Tc ~ 0.1 ms and Tg = 0.1 s. 



5.2. Fixed-potential ICEO 

In the above examples, we have assumed a conducting element which is electrically iso- 
lated from the driving electrodes, which constrains the total charge on the conductor. 
Another possibility involves fixing the potential of the conducting surface with respect to 
the driving electrodes, which requires charge to flow onto and off of the conductor. This 
ability to directly control the induced charge allows for a wide variety of microfluidic 
pumping strategies exploiting ICEO. 
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Material properties of aqueous solution 

Dielectric constant ^ 80eo 7 x 10^^ g cm/V^s^ 

Viscosity rj 10~^ g/cm s 

Ion Diffusivity A 10"^ cm^/s 

Experimental Parameters 

Screening Length A_d 10 nm 

Cylinder radius a 10 /um 

Applied field Eo 100 V/cm 

Electrode Separation L 1 cm 

Characteristic Scales 

Slip velocity Uo = EwEQa/i] 0.7 mm/s 

Cylinder charging time Tc — Xoa/Di 10"'* s 

Electrode charging time Te = XoL/Di 10"* s 

Dimensionless Surface Potential ^ = eEoa/ksT 3.9 

Dukhin number Du 10"^ 



Table 2. Representative values for induced-charge electro-osmosis in a microfiuidic device. 



Perhaps the simplest example of fixed-potential ICEO involves a conducting cylinder 
of radius a which is held in place at a distance h from the nearest electrode (figure 
For simplicity, we consider a ^ h and h <^ L. We take the cylinder to be held at some 
potential Vc, the nearest electrode to be held at Vq, the other electrode at \^ = 0, and 
assume the electrode charging time Tg to be long. In this case, the bulk field (unperturbed 
by the conducting cylinder) is given simply hy Eq = Vq/ L, and the 'background' potential 
at the cylinder location is given by 0(/i) — Vo{l — h/ L). In order to maintain a potential 
Vc, an average zeta potential, 

a ^V,-vJl-y], (5.4) 



is induced via a net transfer of charge per unit length of Xi ~ 27r£mKaCj, along with an 
equally and oppositely charged screening cloud. 

This induced screening cloud is driven by the tangential electric field in the 

standard way, giving a fixed-potential ICEO flow with slip velocity 

FP , rt^'W 



= u. + 2^^ sine [V, -Vo + ^j 0, (5.5) 

where Ug is the quadrupolar 'fixed-total-chargc' ICEO flow Us from (|3.7() and H3.8() . with 
Eo = Vo/L. Note that both the magnitude and direction of the flow can be controlled by 
changing the position h or the potential Vc of the inducing conductor. A freely-suspended 
cylinder would move with an electrophoretic velocity 

away from the position he ~ i(l— K/Vb) where its potential is equal to the (unperturbed) 
background potential. The velocity scale is the same as for fixed-total-charge ICEO, 
although the cylinder-electrode separation h, rather than the cylinder radius a, provides 
the geometric length scale. Since typically h ^ a, fixed-potential ICEO velocities are 
larger than fixed-total-charge ICEO for the same field. 

To hold the cylinder in place against Ue, however, a force is required. Following 
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Figure 4. Fixed-potential ICEO. a) A cylinder of radius a is held a distance h ^ a from a 
nearby electrode and held at the same electrostatic potential Vo as the electrode. A second elec- 
trode is located a distance L away and held at zero potential, so that a field iJo ~ Vb/i is estab- 
lished, b) Leading-order fixed-potential ICEO flow for the system on the left. Like fixed-charge 
ICEO, a non-zero steady flow can be driven with an AC voltage. 



Jeffrey fc Onishi (1981)| the force Fe is given approximately by 



^ \og[ih + Vh^ -a^)/a] -Vh^-aya ^ log (2 ^a) - 1 ' 

and is directed toward he- The fixed-potential ICEO flow around a cylinder held in place 
at the same potential as the nearest electrode (Vc = Vq, he = 0) is shown in figure 01 
The leading-order flow consists of a Stokeslet of strength Fe plus its images, following 



Liron & Blake (1981) The (quadrupolar) fixed-total-charge ICEO flow exists in addition 
to the flow shown, but is smaller by a factor a/h and is thus not drawn. 

With the ability to actively control the potential of the 'inducing' surface, fixed- 
potential ICEO flows afford significant additional flexibility over their fixed-total-charge 
(and also colloidal) counterparts. Note that in a sense, there is little distinction be- 
tween 'inducing' conductors and blocking electrodes. Both impose voltages, undergo 
time-dependent screening, and may drive ICEO flows. Furthermore, their sensitivity to 
device geometry and nonlinear dependence on applied fields open intriguing microfluidic 
possibilities for ICEO flows - both fixed-total-charge and fixed-potential. 



5.3. Simple microfluidic devices exploiting ICEO 

Owing to the rich variety of their associated phenomena, ICEO flows have the potential 
to add a significant new technique to the microfluidic toolbox. Below, we present several 
ideas for microfluidic pumps and mixers based on simple ICEO flows around conducting 
cylinders. The devices typically consist of strategically placed metal wires and electrodes. 
As such, they can be easily fabricated and operate with no moving parts under AC applied 
electric fields. AC fields have several advantages over DC fields: (i) electrode reactions 
are not required to apply AC fields, thus eliminating the concentration and pH gradients, 
bubble formation and metal ion injection that can occur with DC fields, and (ii) strong 
fields can be created by applying small voltages over small distances. For ease of analysis, 
we assume the cylinders to be long enough that the flow is effectively two-dimensional. 
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Figure 5. ICEO micropump designs for 'cross', 'elbow', and 'T' junctions. A conducting cylin- 
der placed in a junction of microchannels, subject to an applied AC or DC field, drives an ICEO 
flow which draws fluid in along the field axis and expels it perpendicular to the field axis. For 
the four-electrode configurations, the field axis can be switched from vertical to horizontal by 
switching the polarities of two diagonally opposite electrodes, reversing the pumping direction. 



5.3.1. Junction pumps 

The simplest ICEO-based devices follow natually from the symmetry of ICEO flows, 
which generally draw fluid in along the field axis and eject it radially. This symmetry 
can be exploited to drive fluid around a corner where two, three, or four micro-channels 
converge at right angles, as shown in figure These DC or AC electro-osmotic pumps 
are reversible: changing the polarity of the four electrodes (shown for the elbow and cross 
junctions in the figure) so as to change the field direction by 90"^, reverses the sense of 
pumping. 

For variety, we give an alternate design for the T-junction pump, which uses a conduct- 
ing plate embedded in the channel wall between the electrodes and cannot be reversed. A 
reversible T junction, similar to the cross and elbow junctions, could easily be designed. 
The former design, however, can be modified to reduce the detrimental effects of viscous 
drag by having the metal plate wrap around to the top and bottom walls (not shown) , in 
addition to the side wall. In general, placing the 'inducing conductor' driving ICEO flow 
on a channel wall is advantageous because it eliminates an inactive surface that would 
otherwise contribute to viscous drag. 

The pressure drop generated by an ICEO junction pump can be estimated on dimen- 
sional grounds. The natural pressure scale for ICEO is P yUa/a ~ EwEq, and the pres- 
sure decays with distance like (a/r)'^. Thus a device driven by a cylinder of radius a in a 
junction with channel half-width W creates a pressure head of order AP ewE^a^ /W'^ . 
For the specifications listed in Table [3 this corresponds to a pressure head on the order 
of mPa. This rather small value suggests that straightforward ICEO pumps are better 
suited for local fiuid control than for driving fluids over significant distances. 



5.3.2. ICEO micro-mixers 



As discussed above, rapid mixing in microfluidic devices is not trivial, since iner- 
tial effects arc negligible and mixing can only occur by diffusion. Chaotic advection 
( Aref (1984) ) provides a promising strategy for mixing in Stokes flows, and various 

Liu et al. (2000)1 



techniques for creating chaotic streamlines have been introduced (e. 



Stroock et al. (2002) ). ICEO flows provide a simple method to create micro- vortices, and 
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Figure 6. AC electro-osmotic mixers. Diffusive mixing in a background flow is enhanced by tfie 
ICEO convection rolls produced by (a) an array of conducting posts in a transverse AC field, 
and (b) conducting objects (or coatings) embedded in channel walls between micro-electrodes 
applying fields along the flow direction. 

L 



h 




Figure 7. Fixed-potential ICEO pumps (left) and mixers (right) can be constructed by elec- 
trically coupling the 'inducing conductor', which may be a post (as shown) or a surface pattern 
(not shown), to one set of electrodes. These devices generalize 'AC electro-osmosis' in flat-surface 
electrodes arrays to other situations, with different flow patterns and frequency responses. 



could therefore be used in a pulsed fashion to create unsteady two-dimensional flows with 
chaotic trajectories. 

In figure Et^, we present a design for an ICEO mixer in which a background flow 
passes through an array of transverse conducting posts. An AC field in the appropriate 
frequency range {t^^ < < is applied perpendicular to the posts and to the mean 
flow direction, which generates an array of persistent ICEO convection rolls. Note that 
the radius of the posts in the mixer should be smaller than shown (a <C W) to validate 
the simple approximations made above, where a 'background' field is applied to each 
post in isolation. However, larger posts as shown could have useful consequences, as the 
final field is amplified by focusing into a smaller region. We leave a careful analysis of 
such issues for future work. 

The same kind of convective mixing could also be produced by a different design, 
illustrated in figure EId, in which an AC (or DC) field is applied along the channel with 
metal stripes embedded in the channel walls. As with the posts described above, the 
metal stripes are isolated from the electrodes applying the driving field. This design has 
the advantage that it drives fiow immediately adjacent to the wall, which reduces 'dead 
space'. 

5.3.3. Devices exploiting fixed-potential ICEO 

By coupling the potential of the posts or plates in the devices above to the electrodes, 
fixed-potential ICEO can be exploited to generate net pumping past the posts. For ex- 
ample, in the cross-junction pump of figurcEl if the central post were grounded to one of 
the pairs of electrodes, there would be an enhanced flow sucking fluid in from the channel 
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between the electrode pair (with an AC or DC voltage). Likewise, a fixed-potential ICEO 
linear pump can be created in the middle of a channel, as shown in figure UK- 

To estimate the flow generated by the single-post device, we note that the post in 
figure if freely suspended, would move with velocity Ue ~ SwVqH/tjL'^ down the 
channel. The post is held in place, however, which requires a force (per unit length) ~ 
47r?7/7£;/(log(W^/a) — 0.9) ( Happel fc Brenner (1983) \. The resulting flow rate depends on 
the length of the channel; however, the pressure required to stop the flow does not. This 
can be estimated using a two-dimensional analog of the calculation of Brenner (1958) for 
the pressure drop due to a small particle in a cylindrical tube, giving a pressure drop 
AP ~ 3F||/4VF. This pressure drop is larger than that for the fixcd-total-charge ICEO 
junction pumps by a factor of O [hW^ / L^) . Furthermore, multiple fixed-potential 
ICEO pumps could be placed in series. When the post is small (a ^ W,L), these pumps 
operate in the same frequency range < w < r~^, as the junction pumps above. 

Clearly, many other designs are possible, which could provide detailed flow control for 
pumping or micro- vortex generation. For example, fixed-potential ICEO can also be used 
in a micromixer design, as shown in figure[33, wherein rolls the size of the channel can be 
established. An interesting point is that frequency response of the ICEO flow is sensitive 
to both the geometry and the electrical couplings. We leave the design, optimization, and 
application of real devices for future work, both experimental and theoretical. 

We close this section by comparing these new kinds of devices with previous examples 
of ICEO in microfluidic devices, which involve quasi-planar micro-electrode arrays to 
produce 'AC electro-osmosis' (Ramos et al. 1999; Ajdari 2000). As noted above, when 
the potential of the 'inducing conductor' (post, stripe, etc.) is coupled to the external 
circuit, it effectively behaves like an electrode, so flxcd-potential ICEO is closely related 
to AC electro-osmosis. Of course, it shares the same fundamental physical mechanism, 
which we call 'ICEO', as do related effects in polarizable colloids that do not involve 
electrodes. Fixed-potential ICEO devices, however, represent a significant generalization 
of AC electro-osmotic planar arrays, because a nontrivial distinction arises between 'elec- 
trode' (applying the field) and 'inducing conductor' (producing the primary ICEO fiow) 
in multi-dimensional geometries. This allows a considerable variety of fiow patterns and 
frequency responses. In contrast, existing devices using AC electro-osmosis peak at a 
single frequency and produce very similar flows (Ramos et al. 1998, 1999; Brown et al. 
2001; Studer et al. 2002; Mpholo et al. 2003). 



6. Surface contamination by a dielectric coating 

The above examples have focused on an idealized situation with a clean metal surface. 
In this section, we examine the effect of a non-conducting dielectric layer which coats the 
conductor, and find that any dielectric layer which is thicker than the screening length Au 
significantly reduces the strength of the ICEO flow. Furthermore, the ICEO flow around 
a dielectric object, rather than a perfectly conducting object as we have discussed so far, 
is presented as a limiting case of the analysis in this section. 

We start with a simple physical picture to demonstrate the basic effect of a thin 
dielectric layer. Consider a conducting cylinder of radius a coated with a dielectric layer 
of thickness Xd ^ a (so that the surface looks locally planar) and permeability e^, as 
shown in flgure|Sl In steady state, the potential drop from the conducting surface </) = 
to the potential 0oo outside of the double layer occurs across in two steps: across the 
dielectric (where E = Ed), and across the screening cloud (where E = E^,), so that 



EdXd + EyjXo — 0OO. 



(6.1) 
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Figure 8. A dielectric layer of thickness Ad and permittivity Sd coating a conductor. The 
potential drop between the external potential <j)ao and at the conductor occurs in two steps: 
A(l)d ~ Ed\d across the dielectric and Afjiu, ~ Ew^d across the double layer. 



The electric fields in the double layer and in the dielectric layer are related via SdEd 



EyjEyj, so that 



1 + — — ^ ) EyjXo = (PoD- 
£d Ad, 



(6.2) 



Since E^iXd is approximately the potential drop ( across the double-layer, and since 
the steady-state bulk potential is given by H3.4|l to be (/)oo = 2EQacos9, we find the 
induced-charge zeta potential to be 



c 



2Eoa cos 9 

1 + Sw^d/Sd^D 



(6.3) 



Thus unless the layer thickness Xd is much less than edXa/sw, the bulk of the potential 
drop (jjoo occurs across the dielectric layer, instead of the double-layer, resulting in a 
reduced electro-osmotic slip velocity. 

The modification to the charging time Tc for the coated cylinder can likewise be un- 
derstood from this picture. The dielectric layer represents an additional (parallel-plate) 
capacitor of separation Xd and filled with a dielectric Ed, in series with the capacitive 
screening cloud, giving a total capacitance 



and a modified RC time 



Apa 



1 



rl An 



£d^D 



£'il! Xn 



£dA 



D 



(6.4) 



(6.5) 



as found by Ajdari (2000) for a similar calculation for a compact (Stern) layer. A dis- 
cussion of double-layer capacitance in the nonlinear regime, where the compact layer is 



approximated by a thin dielectric layer is given by Macdonald (1954) 



A full calculation of the induced zeta potential around a dielectric cylinder of radius a 
coated by another dielectric layer with outer radius b and thickness (b — a) is straightfor- 
ward although cumbersome. (Note that the analogous problem for a coated conducting 
sphere was treated by A. S. Dukhin (1986) as a model of a dead biological cell.) The 



resulting induced-charge zeta potential is 



c = 



2&£'o(l + rc)cos6' 
2 + Kb{l-r^) ' 



(6.6) 
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where Tr is defined to be 



62 + a2 + £,^/erf(62 „ ■ 
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(6.7) 



(6.8) 



It is instructive to examine limiting cases of the induced zeta potential (|6.6() . In the 
limit of a 'conducting' dielectric coating ed/£w oo, we recover the standard result 
for ICEO around a metal cylinder, as expected. In the limit of a thin dielectric layer 
6 = a + Xd, where ^ a, the induced zeta potential is given by 

2bEQ cos 9 



C(Ad < a) 
as found in H6.3|l . with a charging time 

Tc{Xd < a) w 



1 + Xd£w/ XoEd ' 



Xpb 
D 



£wXd 

SdXp 



(6.9) 



(6.10) 



as found in H6.5|l . Therefore, the ICEO slip velocity around a coated cylinder is close 
to that of a clean conducting cylinder (C — 2bEo cos 9) only when the dielectric layer is 
much thinner than the screening length times the dielectric contrast, Xd <C XoSd/sw 
The zeta potential induced around a conducting cylinder with a thicker dielectric layer, 



({Xd > Xd) 



2bEo — ^ cos ( 



(6.11) 



(6.12) 



£wXd 

is smaller by a factor of O (Ad/A^), and the charging time, 

M \ \ ^ £wXd Xob 

Tc{Xd > Xd) ~ — ^ p-, 

£dXd D 

is likewise shorter by a factor of O (X^/Xd)- As a result, strong ICEO flow requires a 
rather clean and highly polarizablc surface, with minimal non-conducting deposits. 

Note that the limit of a pure dielectric cylinder of radius b is found by taking the limit 
a — !■ 0. This gives a zeta potential of 



({a = 0) 



2bEo cos 9 



1 



2—EqXd cos 9, 



(6.13) 



U b/edXo 

and an ICEO slip velocity which is smaller than the conducting case (|3.7f) by O (Ad/6). 
The charging time for a dielectric cylinder is given by 



P X2 

Tc{a = 0) « — , 

£d D 



(6.14) 



as expected from H6.5(> in the limit A^ ^ Xd- 



7. Induced-charge electro-osmosis: systematic derivation 

In this section, we provide a systematic derivation of induced-charge electro-osmosis, in 
order to complement the physical arguments above. We derive a set of effective equations 
for the time-dependent ICEO flow around an arbitrarily-shaped conducting object, and 
indicate the conditions under which the approximations made in this article are valid. 
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Starting with the usual 'electrokinetic equations' for the electrostatic, fluid, and ion fields 
(as given by, e.g. , Hunter (2000)| ), we propose an asymptotic expansion that matches an 



inner solution (valid within the charge cloud) with an outer region (outside of the charge 
cloud), and that accounts for two separate time scales-the time for the charge cloud to 
locally equilibrate and the time scale over which the external electric field changes. 

Although it has mainly been applied to steady-state problems, the method of matched 
asymptotic expansions is well established in this setting, where it is commonly called 
the 'thin double layer approximation'. For simplicity, like many authors, we perform our 
analysis for the case of a symmetric binary electrolyte with a single ionic diffusivity in a 
weak applied field, and we compute only the leading-order uniformly valid approximation. 



The same simplifying assumptions were also made by Gonzalez et al. (2000) in their 
recent asymptotic analysis of AC electro-osmosis. 

In the phenomenological theory of diffuse charge in dilute electrochemical systems, the 
electrostatic field obeys Poisson's equation, 



V20=-(^^i^^, (7.1) 

where n± represent the local number densities of positive and negative ions. These obey 
conservation equations 

^ + V.(n±v±)=0, (7.2) 

where v± represent the velocities of the two ion species, 

v± = =F&eV0 - keThV log n± + u, (7.3) 

where b is the mobility of the ions in the solvent and u is the local fluid velocity. Terms 
in H7.3|l represent ion motion due to (a) electrostatic forcing, (b) diffusion down density 
gradients, and (c) advcction with the local fluid velocity. The fluid flow obeys the Stokes 
equations, with body force given by the product of the charge density with the electric 
field, 

r^V^u - Vp = e(n+ - ) V0, (7.4) 

along with incompressibility. Strictly speaking, this form does not explicitly include os- 
motic pressure gradients, as detailed below. However, osmotic forces can be absorbed 
into a modified pressure field p, and the resulting flow is unaffected. 

For boundary conditions on the surface F of the conductor, we require the fluid flow to 
obey the no-slip condition, the electric potential to be an equipotential (with equilibrium 
zeta potential ^o), and the ions to obey a no-flux condition: 

u(r) = (7.5) 
0(F) = Co (7.6) 
fi-v±(F) = 0, (7.7) 

where fi represents the (outer) normal to the surface F. Far from the object, we require 
the fluid flow to decay to zero, the electric fleld to approach the externally-applied electric 
field, and the ion densities n± to approach their constant (bulk) value np. 

In order to simplify these equations, we insert (|7.3(l in (|7.2(l . and take the sum and 
difference of the resulting equations for the two ion species to obtain 

cP + Dk^cP - DV'^cP - ehV ■ [c^Vc/)] + u • Vc^ = (7.8) 
c" - DV^c" - ebV ■ [cPV<j)] + u • Vc" = 0, (7.9) 
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where we have used H7.1|l and where we have defined 

c"" = n+ + n-- 2no (7.10) 
cP^n+-n^. (7.11) 

The first variable, c'^ , represents the excess total concentration of ions, while the second, 
qP, is related to the charge density via c'' = p/e. The first three terms in (|7.8|) represent 
a (possible) transient, electrostatic transport, and diffusive transport, respectively, and 
will be seen to give the dominant balance in the double layer. The next term represents 
the divergence of a flux of excess ionic concentration c'^ (but not excess charge) due 
to an electric field. The final term represents ion advection with the fluid flow u. The 
analogous H7.9() for the charge-neutral c^ lacks the electrostatic transport term. Both c*^ 
and c'' decay away from the solid surface and obey the no-flux boundary conditions (|7.7|) , 

fi-Vc''|r = -r^(2no + c^)n- V0|r (7.12) 
fi- Vc^lr^-T^Cfi-V^lr, (7.13) 

at the surface F. 

In what follows, we obtain an approximate solution to the governing equations (|7.1|) . 
H7.8|) . and (|7.9|) . at the leading order in a matched asymptotic expansion. We first analyze 
the solution in the 'inner' region within a distance of order Xd of the surface. Non- 
dimensionalization yields a set of approximate equations for the inner region and a set of 
matching boundary conditions which depend on the 'outer' solution, valid in the quasi- 
neutral bulk region (farther than Xo from the surface). Similarly, approximate equations 
and effective boundary conditions for the outer region are found. By solving the inner 
problem and matching to the outer problem, a set of effective equations are derived for 
time-dependent, bulk ICEO flows. 

7.1. Non-dimensionalization and 'inner' region 

We begin by examining the 'inner region', adjacent to the conducting surface. Denoting 
non-dimensional variables in the inner region with tildes, we scale variables as follows, 

r = Xdv, t = {k^D)-H, (j) = $00, u = UqU, c^ = 2no^dP, c" = 2?^o^'^c^ E = EqE 

(7.14) 

where $0 a-nd Uq are potential and velocity scales (left unspecified for now), and where 
we have introduced the dimcnsionless surface potential. 

Note that c*^ scales with ^'^ to satisfy the dominant balance in H7.9|) . Finally, we have 
scaled time with the Debye time, to = {k?D)~^ = X\)/D, although the analysis will 
dictate another time scale, as expected from the physical arguments above. 
The dimcnsionless ion conservation equations H7.8I7.9|I become 

dcP 



dt 



dt 

where we have introduced the Peclet number. 



c^Wct) 
cPV^ 



-Peu-Vc'' = (7.16) 
-Peu-Vc'==0, (7.17) 



Pe^i|. (7.18) 
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The boundary conditions (|7.12I7.13|I in non-dimensional form are given by 

n- VS^Ir ^-Cl + ^-^e")!!- V0|r (7.19) 
n- Vc^lr = -c''fi- V0|r. (7.20) 

In the analysis that follows, we concentrate on the simplest limiting case. We assume 
the screening length to be much smaller than any length Lq associated with the surface 
geometry, parametrized through 

e - (AiLo)-' = Ap/Lo « 1, (7.21) 

so that the screening cloud 'looks' locally planar. As mentioned above, the singular 
limit of thin double layers, e ^ 1, is the usual basis for the matched asymptotic ex- 
pansion. The regular limit of small Peclet number, Pe ^ 1, which holds in almost any 
situation, is easily taken by setting Pe = 0. Finally, we consider the regular limit of 
small (dimensionless) surface potential, ^' <C 1, which is the same limit that allows the 
Poisson-Boltzmann equation to be linearized. With these approximations, the system is 
significantly simplified: First, cP is coupled to only through terms of order O {^^) in 
H7.16|l and boundary condition H7.19|l . Second, is smaller than cP by a factor 5*, and 
is thus neglected: = {) + O (*^). 

In this limit, we obtain the linear Debye equation and boundary condition for cf alone: 

^ + 5" = e2y2~p (7 22) 
dt 

n- V£''|r = -fi- V^lr. (7.23) 
The potential is then recovered from (|7.1|l in the form, 

V^^ = -cP, (7.24) 

with the far-field boundary condition, 

and the (equipotential) surface boundary condition 

0(r) = $oCo. (7.26) 

where C,q is the dimensionless equilibrium zeta potential. 

Since e ^ 1, we introduce a locally Cartesian coordinate system {n, Z}, where n is 
locally normal to the surface, and I is locally tangent to the surface. The governing 
equations H7.22|l and (|7.24|) are both linear in and 0, which allows the electrostatic 
and ion fields to be expressed as a simple superposition of the equilibrium fields, 

= -<^cq - Coe-", (7.27) 

and the time-dependent induced fields c^* and (pi, via 

c'^~cP^ + c^ (7.28) 

= 0cq + 0i. (7.29) 

The induced electrostatic and ion fields must then obey 

72 



V"./), = -cf (7.30) 

dt 



^ + cf-V25f = 0, (7.31) 
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subject to boundary conditions 



^,{T) = (7.32) 
V5f|r = -ft-V0,|r, (7.33) 



on the surface. 

Another consequence of the Unearity of (|7.22|l and (|7.24|) . and therefore of the decou- 
phng of the equihbrium and induced fields, is that different characteristic scales can be 
taken for the two sets of fields. We scale the potential in the equilibrium problem with 
the equilibrium zeta potential ^q, and the potential in the induced problem with EqLq, 
the potential drop across the length scale of the object. In order that the total surface 
potential be small (^ ^ 1), however, we require that the total zeta potential be small, 

M±l^ « 1. (7.34) 

At the end of this section, we will briefly discuss the rich variety of nonlinear effects 
which generally occur, in addition to ICEO, when this condition is violated. 

We expect the charge cloud c^ and electric potential (p to vary quickly with n, but 
slowly with I (along the surface). Furthermore, we expect the charge cloud to exhibit two 
time scales: a fast (transient) time scale, over which reaches a quasi-steady dominant 
balance, and a slow time scale Tc, over which the quasi-steady solution changes. The 
physical arguments leading to H4.6|l for the charging time suggest that this slow time 
scale is given by Tc ~ 1/e, which is not obvious a priori, but which will be confirmed by 
the successful asymptotic matching. 

In order to focus on the long-time dynamics of the induced charge cloud, and guided 
by the above expectations, we attempt a quasi-steady solution to (|7.3()|l - H7.33(l of the 
form 

cf = cf(n,ef,et) (7.35) 
4>i = 4>iin, el, ei). (7.36) 

It can be verified that 

cf ^ A{el ei>-" - e^li^Lf^lne-" (7.37) 

4>^ = -Aiel, et)e-" + [26^" + ne""] + B{el, et) + C{el, et)n, (7.38) 

solve the governing equations to O (e^). 

The equipotential boundary condition (|7.32|l is satisfied when 

A{el d) = S(e[, ei) ^ eA{el, ei), (7.39) 

so that the induced double-layer charge density is proportional (to leading order in e) to 
the potential just outside the double layer. The normal ion flux condition l|7.33(l . 

eA{elei) ^ C{d,ei), (7.40) 

relates the evolution of the double-layer charge density to the electric field normal to the 
double-layer. Matching the inner region to the outer region provides the final relations. 
In the limit n — > oo, the electrostatic and ion fields approach their limiting behavior 

cP ^0 (7.41) 
4> B{el, ei) + C{el, et)fi. (7.42) 
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7.2. Outer solution 

We now turn to the region outside of the double-layer, and use overbars to denote 'outer' 
non-dimensional variables (scaled with length scale Lq and potential scale EqLq). Ac- 
cording to H7.37|) and (|7.41|) . any non-zero charge density that exists within the inner 
region decays exponentially away from the surface T. A homogeneous solution {c^ = 0) 
thus satisfies (I7.22|) and the decaying boundary condition at infinity. With eg = 0, H7.24|l 
for the electrostatic field (/)o in the outer region reduces to Laplace's equation, 

V'^o = 0, (7.43) 

with far-field boundary condition (|7.25|l given by 

V0O ^ -E. (7.44) 

To determine (j)o uniquely, one further boundary condition on the surface T is required, 
which is obtained by matching to the inner solution. The limiting value of the 'outer' 
field <j)o, 

^o(r^r)^0o(r) + ^±(r)n, (7.45) 

must match the inner solution (|7.42|) . which gives two relations 

Biel,ei) = Ml), (7.46) 
C{el,ei) ^eE^{l). (7.47) 

Using l|7.47(l in l|7.4()|l . the fact that C is C (e) verifies that the assumed time scale 
Tc = Tu/e for induced double layer evolution is indeed correct, as expected from physical 
arguments (|4.6() . 

7.3. Effective equations for ICEO in weak applied fields 

Using the above results, we present a set of effective equations for timc-dcpcndcnt ICEO 
that allows the study of the large-scale flows, without requiring the detained inner so- 
lution. What emerges is a first-order ODE for the dimensionless total surface charge 
density in the diffuse double layer, q = A, so an 'initial' value for q must be specified. 
From H7.39|l and H7.46|l . the 'outer' potential on the surface V is given by 

0o(r,ei)=g(r,e<), (7.48) 

which, along with the far-field boundary condition l|7.25(l . uniquely specifies the solution 
to Laplace's equation l|7.43(l . From this solution, the normal field i?j_(r, et) is found, 
which (using H7.47|l and (j7.40|l ^ results in a time-dependent boundary condition, 

|| = 6^^(r,ei), (7.49) 
which is more naturally expressed as 

|| = ^±(r,t), (7.50) 

in terms of the dimensionless time variable, t = et = t/Tc. 

We have thus matched the bulk field outside the charge cloud with the 'inner' behavior 
of the charge cloud in a self-consistent manner. The 'inner' solutions for and C, equi- 
librate quickly in response to the (slow) charging, and affect the boundary conditions 
which determine the 'outer' solution. Matching (|7.47() and (|7.40() results in a relation 
between the charge cloud and normal ionic flux, confirming the validity of 1)3. 3|l and 
H4.3|l . which were previously argued in an intuitive, physical manner. This analysis has 
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demonstrated that errors to this approach are of order O (^^), O (P) and O (e^)- In 
this Umit, the system evolves with a single characteristic time scale, Tc = XdLo/D, set 
by asymptotic matching, which physically corresponds to an RC coupling, as explained 
above. 

7.4. Fluid dynamics 

7.4.1. Fluid body forces: electrostatic and osmotic 

Thus far, we have derived the effective equations for the dynamics of the induced 
double-layer. To conclude, we examine the ICEO slip velocity that results from the 
interaction of the applied field and the induced diffuse layer. We demonstrate that, in the 
limits of thin double layer and small surface potential, the Smoluchowski formula l|2.4|l 
for fluid slip holds. 

Following Levich (1962) (p. 484), we consider a conducting surface F immersed in a 



fluid, with an applied 'background' electrostatic potential 0o. In addition, a thin double- 
layer (with potential ^) exists and obeys 

^ _ (n+-n_)e ^ _noe _ ^^^^^^ 

so that the total potential (f) = (t>o + Note that this assumes that the double layer is in 
quasi-equilibrium. 

Fluid stresses have two sources: electric and osmotic. The electric stress in the fluid is 
given by the Maxwell Stress tensor, 

= e„, [e,E, - • , (7.52) 

from which straightforward manipulations yield the electrical body force on the fluid to 
be 

Fb = V • T = £^V(/)V20 = e^V(0o + (7.53) 
Osmotic stresses come from gradients in ion concentration, and exert a fluid body force 

Fo = -keTV [n+ + n_) = eVC {n+ - n_) = -e^VS,V^t (7.54) 

Thus the total body force on the fluid, is given by the sum of F^; and Fq, 

F = £.„V(/)oV2e = pEB. (7.55) 

Therefore, when both electric and osmotic stresses are included, the body force on the 
double layer above a conductor is given by the product of the local charge density p and 
the 'background' electric field E^ = — V</>o (which, importantly, docs not vary across the 
double layer). Note also that the same fluid flow would result if the double-layer forcing 
p(EiB + E?) were to be used, since the osmotic component ()7.54() is irrotational and can 
be absorbed in a modified fiuid pressure. 

7.4.2. ICEO slip velocity 

Finally, we examine the ICEO slip velocity that results when the electric field — V(/)o 
drives the ions in the induced charge screening cloud p. 

We look first at the flow in the 'inner' region of size Ajj. For the ICEO slip velocity to 
reach steady state, vorticity must diffuse across the double layer, which requires a very 
small time t^{\d) = ^d/^ ~ 10" s. Because t^^ is so much faster than the charging 
time and Debye time, we consider the ICEO slip velocity to follow changes in ^ or ^ 
instantaneously. As noted above, the unsteady term may play a role in cutting off the 
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oscillating component of the bulk flow. However, our main concern is with the steady, 
time-averaged component, which simply obeys the steady, forced Stokes equation. 

Non-dimensionalizing as above, we re-express the forced Stokes equations (|7.4|) . with 
forcing given by H7.55|l . using a stream function ij} defined so that ui = dn4' ^^nd u„ = 
—diip. The stream function obeys 

VV^^fe^-S^V (7.56) 

\dn dl dl dnj ' ^ ^ 

where the stream function has been scaled by t/j = (eiuE'oCo/^k)?/'. We perform a local 
analysis around the point Z = of the ICEO flow driven by an applied tangential electric 
field {-Ell, E±}. Using representations of and (j) around I = 0, 

cP = cP{el)e-^ (7.57) 
^ = -E\\l{el, en) - E^n{el, en), (7.58) 

it is straightforward to solve H7.56|l to O (e). The tangential and normal flows are then 
given by 



Ul = 



(^Cm\ (0 + e^E^^ (1 - e-™) + (3 - (3 + ^n)e-™) 



(7.59) 



-e- 



-(Ci?||)(l-«n-e-™), (7.60) 
M ol 

where C(0^|| contains terms of O (e). To leading order, then, the slip flows obey 

ui^ ~—C{l)Ei\{l)+0{e) and w„ = 0(e). (7.61) 

The tangential flow does indeed asymptote to H2.4|l . with local zcta potentials and back- 
ground field, and the normal flow velocity is smaller by a factor of order O (e). 

Thus an ICEO slip velocity is very rapidly established in response to an induced zeta 
potential and 'outer' tangential field. Furthermore, despite double-layer and tangential 
field gradients, the classical Helmholtz-Smoluchowski formula, (|2.4|l . correctly gives the 
electro-osmotic slip velocity. This may seem surprising, given that the tangential field 
vanishes at the conducting surface. 

The final step involves finding the bulk ICEO flow must be found by solving the 
unsteady Stokes equations, with no forcing, but with a specified ICEO slip velocity on 
the boundary T, given by solving the effective clcctrokinetic transport problem above. 

7.5. Other nonlinear phenomena at large voltages 

Although we have performed our analysis in the linearized limit of small potentials, it 
can be generalized to the 'weakly nonlinear' limit of thin double layers, where H2.7|) holds 
and Du <C 1. In that case, the bulk concentration remains uniform at leading order, and 
the Helmholtz-Smoluchowski slip formula remains valid. The main difference involves the 
time-dependence of double-layer relaxation, which is slowed down by nonlinear screening 
once the thermal voltage is exceeded. It can be shown that the linear time-dependent 
boundary condition, (|3.3(l or H7.50|l . must be modified to take into account the nonlinear 
differential capacitance, as in H4.12|l . Faradaic surface reactions and the capacitance of 
the compact Stern layer may also be included in such an approach, as in the recent work 
of [Bonnefont et al. (200T]1 



The condition l|^.7|l tor the breakdown of Smoluchowski's theory of electrophoresis 
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(with C = Co + Equ) coincides with the condition Tc(C) <C td, where Tc takes into ac- 
count the nonhnear differential capacitance H4.14|l . When this condition is violated (the 
'strongly nonlinear' regime), double-layer charging is slowed down so much by nonlinear- 
ity that it continues to occur at the time scale of bulk diffusion. At such large voltages, 
the initial charging process draws so much neutral concentration into the double layer 
that it creates a transient diffusion layer which must relax into the bulk, while coupled 
to the ongoing double-layer charging process. Bazant, Thornton & Ajdari (2004) explore 
such processes. 

In the strongly nonlinear regime, the Dukhin number is typically not negligible, and 
bulk concentration gradients (and their associated electrokinetic effects) are produced by 
surface conduction (see S. S. Dukhin (1993) and Lyklema (1991) I. Tangential concentra- 
tion gradients modify the usual electro-osmotic slip by changing the bulk electric field 
(concentration polarization) and by producing diffusio-osmotic slip. Therefore, both the 
steady state and the relaxation processes for ICEO flows arc affected. 

Finally, large voltages can also lead to the breakdown of ideal polarizability via sponta- 



neous Faradaic reactions on different sides of the object, as noted by Murtsovkin (1991) 



Gamayunov et al. (1992) observed that the induced-charge electro-osmotic flow around 
metal colloidal spheres reverses direction for large colloids, and argued that Faradaic 
reactions were responsible. Furthermore, Barany et al. (1998) measured large conduct- 
ing colloids to have 'superfast' (second-kind) electrophoretic velocities. They argued that 
sufficiently strong fields cause Faradaic currents at the two sides of the colloid (as though 
they were electrodes), resulting in the development of a bulk 'space charge' and, corre- 
spondingly, second-kind electrophoresis. 

A complete description of time-dependent ICEO at large voltages, which is beyond the 
scope of this article, would require considering all of these effects at once. The approx- 
imation of thin double layers, which has been applied mostly to steady-state problems 
involving non-polarizable objects, is a good starting point. However, the presence of mul- 
tiple length and time scales complicates mathematical analysis, especially in any attempt 
to go beyond the leading-order approximation. An important goal, therefore, would be to 
extend the method of matched asymptotic expansions to derive effective equations and 
boundary conditions for strongly nonlinear ICEO and to carefully analyze asymptotic 
corrections. 



8. Summary and Discussion 

In this article, we have described the general phenomenon of induced-charge electro- 
osmosis (ICEO), which includes a wide variety of techniques (both old and new) for 
driving steady micro-flows around conducting or dielectric surfaces using AC or DC elec- 
tric fields. Wc have given a physical picture of the basic mechanism for ICEO, involving 
the inhomogcneous surface charge induced in the conductor in order to maintain an 
equipotential surface in the presence of an applied field. In response, the electric field 
normal to the surface/charge cloud drives ions into an inhomogcneous (dipolar) charge 
cloud, which are in turn driven by the tangential electric field. This results in ICEO slip 
velocities of magnitude Uq ~ EwEQa/rj. A charging time scale Tc ^ Xd<i/ D is required for 
these induced charge clouds to form. Due to the dependence on the square of the applied 
field Eq, a nonzero time-averaged ICEO flow can be driven using AC fields of sufficiently 
low frequencies (w <C I/tc). 

We have performed explicit calculations for the steady and unsteady ICEO slip veloci- 
ties (suddenly-applied and sinusoidal fields) around symmetric conducting cylinders. The 
ICEO flow for conducting cylinders is quadrupolar and decays with distance like 7'"^. We 
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have also performed a systematic, matched-asymptotic analysis of the equations for the 
ion transport, electrostatics and fluid flow to confirm the validity of the physically intu- 
itive approach. The analysis produces a set of effective equations which 'integrate out' 
the dynamics of the thin screening cloud, allowing the bulk ICEO flow to be calculated 
with macro-scale calculations alone. 

We have also considered polarizable dielectric surfaces for two reasons. First, we have 
shown that a dielectric layer of thickness Ad can reduce the strength of the ICEO slip 
velocity by a factor of order Xo/a when the dielectric layer is sufficiently thick (Ad ^ A_d). 
This underscores the necessity of using clean and/or treated surfaces to ensure a clean 
conductor/ water interface. Second, an ICEO flow is set up even around a purely dielectric 
colloidal particle of permittivity e^, but with reduced ICEO slip velocity, [/q ^ SdE^Xo M- 

In this article, we have concentrated upon ICEO flows in systems of high symmetry: 
circular cylinders and spheres in spatially uniform applied fields, for which simple exact 
solutions are possible. In a companion article, we will explore the implications of broken 
spatial symmetries - both via asymmetric surface properties and gradients in the applied 
electric field - along with more potential applications to microfluidic devices. For a brief 



summary of our results, the reader is referred to Bazant & Squires (2004) 



In conclusion, ICEO is a rather general and potentially useful phenomenon, capable of 
producing large fluid 'slip' velocities around polarizable surfaces, under AC or DC fields. 
Many variants exist on the basic situations presented in this article. For example, one 
can apply spatially inhomogeneous electric fields, vary the geometry or electrical prop- 
erties of the polarizable surface, apply fixed-potential (or actively-controlled potential) 
ICEO flows, and so on. The directions seem promising to pursue experimentally in real 
microfluidic devices. 

In the presence of bulk concentration gradients produced by surface conduction, Faradaic 
processes, or transient double-layer adsorption, more general electrokinetic phenomena 
may also occur at polarizable surfaces. These effects have been described to varying de- 
grees in the Russian literature on 'non-equilibrium electrosurface phenomena' in colloidal 
systems, especially over the past few decades. It is our hope that this mature subject, 
which includes ICEO as a limiting case, will gain renewed attention from the microfluidics 
community in the coming years. 
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